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ABSTRACT 

We present a novel method for simulation of the interior of large cosmic voids, suitable for study of 
the formation and evolution of objects lying within such regions. Following Birkhoff's theorem, void 
, regions dynamically evolve as universes with cosmological parameters that depend on the underdensity 

' of the void. We derive the values of Om, ^^a, and iJg that describe this evolution. We examine how 

I the growth rate of structure and scale factor in a void differ from the background universe. Together 

04 . with a prescription for the power spectrum of fluctuations, these equations provide the initial conditions 

^ ' for running specialized void simulations. The increased efficiency of such simulations, in comparison 

1^ \ with general-purpose simulations, allows an improvement of upwards of twenty in the mass resolution. 

■ As a sanity check, we run a moderate resolution simulation [N = 128^ particles) and confirm that the 

resulting mass function of void halos is consistent with other theoretical and numerical models. 

Subject headings: cosmology: large-scale structure of the universe — cosmology: theory - methods: 
n-body simulations 
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1. INTRODUCTION 



> , 

Large cosmic voids provide a unique testbed for models of galactic formation and evolution (Peebles 2001). In recent 
years, there has been intense effort to identify and study voids in a number of large surveys such as the IRAS 1.2 Jy Survey 
(El- Ad et al. 1997), the PSCz Redshift Survey (PHonis & Basilakos 2002; Hoyle & Vogeley 2002), the Las Campanas 
Redshift Survey (Arbabi-Bigdoh & Miiller 2002) and the Updated Zwicky Catalog (Hoyle & Vogeley 2002). Observed 
' large voids have typical diameters of D ~ 25/i~^Mpc, density contrast 5 ~ —0.9, and fill at least 40% of the universe 
I (Hoyle & Vogeley 2002). The properties of galaxies in voids clearly differ from those in higher-density regions (Grogin & 
■ Geller 1999a, b; Rojas et al. 2003a, b). Recent studies using samples of 10^ void galaxies from the Sloan Digital Sky Survey 
' demonstrate that void galaxies are bluer, fainter, more disklike, and have higher specific star formation rates (Rojas et 
^T^: al. 2003a,b). 

Q A variety of simulation methods has been applied to study formation of object in voids, including application of 
\ Semi- Analytic Models (SAMs) to compute properties of galaxies associated with dark matter halos in cold dark matter 
CZJ . simulations (Benson et al. 2003; Mathis & White 2002), simulations including hydrodynamics (Ostriker et al. 2003), and 
' specialized high-resolution simulations of dark matter halos in voids (Gottlober et al. 2003). Studying the characteristics 
of galaxies in the context of large-scale numerical simulations can be extremely costly because the region of interest 
represents only ^ 6 — 8% of the total mass of the simulation, and less than half the volume. Thus, enormous amounts 
of computing time must be spent in order to study voids in a particular cosmological model. Studies over many models 
naturally become even more difficult. 

In principle, simulation of voids is far easier than simulation of a typical large patch of the universe, because voids are 
dynamically much simpler systems than clusters. Large voids remain in the linear regime for relatively longer time and 
approach spherical symmetry for much of their evolution. Birkhoff's theorem (Birkhoff 1923) tells us that the internal 
dynamics of a spherically symmetric system will be independent of the dynamics of the outside universe. This theorem has 
been applied in the development of spherically symmetric infall models for galaxies (e.g., Gunn & Gott 1972; Schechter 
1980), cluster infall regions (Regos & Geller 1989), for large-scale inhomogeneities in general (Silk 1974), and for studying 
the population of Lya clouds in voids (Manning 2002). Voids can be approximated as expanding, isolated universes unto 
themselves which do not accrete matter from the universe at large. As we describe below, this is only an approximation, 
but one that is useful for studying formation of objects well within the edges of voids. 

We propose to exploit this simplicity to simulate universes entirely composed of void- like regions. A patch anywhere 
within such a simulation would mimic the interior of a large void and allow investigation of such issues as: What are 
galaxies like in such a universe? At what epoch did they form? 

The organization of this paper is as follows. In §2, we derive the basic relations between the cosmological parameters: 
expansion coefficients, and local values of i^Af, Oa and 77, in the void region and outside of it. In §3, we discuss the 
primordial power spectrum which should be used in running the simulations. Details of applying this approach are 
discussed in §4. In §5, we discuss the results of a small simulation meant as a "sanity check" on the proposed method. 
Finally, in §6, we discuss these results and future applications. 

2. COSMOLOGICAL PARAMETERS IN THE VOIDS 

We can use Birkhoff's theorem to explore the physics of underdense regions, because voids have the properties that 
they are largely isolated throughout their evolution, and they are roughly spherical. Birkhoff's theorem tells us that a 
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spherically symmetric region with specified mean density, cosmological constant density, radius r, and expansion rate f 
will evolve exactly like a universe with a Hubble constant r/r, and a mean density equal to that of the spherical region. 
The growth of perturbations within that region will exactly mimic (up to long-range tidal forces) structure growth in a 
universe with the fiducial cosmology. 

For the purpose of our discussion, we assume a background cosmology with a pure CDM+Baryon matter density of 
Po at the present (throughout, we will use the convention that a subscript "0" denotes the value of the parameter at 
the present). All calculations assume a negligible density of relativistic matter. The matter density parameter may be 
expressed as = SttGpo/BHq. We also assume a known cosmological constant A with corresponding density parameter 
Oa. 

At an arbitrarily high redshift, when the scale factor = a{ti) = 1/(1 + z^), consider a spherical imderdense region 
with physical size r„(fi) (henceforth, subscript "v" denotes a property of the void region). If the density contrast of the 
void is 6v{ti), then the mass contained in the region is, of course, 

M. = ^r,mi + S.it)]^ (1) 

which is a conserved quantity at all times. For reasons of convenience, we also define the dimensionless variable 

which, at very early times, gives a„(f) = a{t). From equations (1,2) we get the relation 
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If the void is spherically symmetric, its evolution is described by the Friedmann equation, 

'ry{t)\'^ 8lTGpy{t) A„ ky 



ry{t)J 3 '3 Ry{tY 



(4) 



where At,, the cosmological constant within the void, is identical to A in the rest of the universe (thus the subscript v will 
be dropped from this term in subsequent discussion). Ry{t) is the (presently unknown) radius of curvature within the 
void, and is the sign of the curvature within the void. 

The right-hand side of equation (4) may be simplified considerably. Equation (1) yields the relation 

Pv{t) = {l + 5^a{t))^ (5) 

From this, we may expand the right-hand side of equation (4) as 



where wc have substituted (a/a)^ for the right-hand side of the Friedmann equation of the background cosmology. Note 
that this substitution assumes the background cosmology is spatially flat. It is a trivial, albeit algebraically messy exercise 
to extend this analysis to the case of non-flat background cosmologies. Note also that in pathologically closed universes, 
the void itself may be closed as well. 

We may also use equation (3) to expand the left-hand side of equation (4) as 
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Prior to equation (7), the relations we derive are fully general. However, the second two equalities in equation (7) only 
hold at extremely early times in the evolution of the void, since at early times dy oc a and \6y \ « 1. 



Equations (6,7) may be combined, yielding (at early times only) 
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Substituting, once again, H'^{t) = Hq{Qm/ci^ + ^a), we obtain 

where we have substituted ky = —1 (necessarily the case for an underdense void in a flat cosmology), or 



It is clear that Ry oc Uy at all times. Thus, by looking at the limiting case of t — > 0, we find at all times 



where the term in the radical is a constant provided ti is selected sufficiently early. This ratio may be given as a fixed 
parameter of the model. As we will see, the ratio 

V = ^ (12) 

together with the parameters of the background cosmology {Qm,^a,Ho) completely constrain the evolution of the void 
region. Wc will compute the final undcrdcnsity of the region, (5„(to), the equivalent density parameter, ^v,M, the equivalent 
cosmological constant, ri„,A, and the equivalent Hubble constant, Hy{t) as a function of the parameter rj. 

Having specified r], we may again write the general form of the Friedmann equation within the void (valid for all times) 

= o 1- IlQilA — -ng - 
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To proceed, wc need to solve a boundary value problem using equation (13). By integrating from <; = to to, and given 
the constraint that at early times, a„(f) = a{t), we solve for the parameter 

a = ay{to) (14) 

Because the void region expands faster than the backgroimd universe, a > 1. Thus, a = a{ri), and wc obtain the 
relationships for the density contrast and effective cosmological parameters in the void at the present epoch, in terms of 
the parameter a: 

Svito) = (15) 



Hy,o = HoJ^ + n^-'4^ (16) 
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where Clk is the corresponding "density" parameter of the curvature term in the Friedmann equation, such that Qm + 

Oa + J^/c = 1- 

In the hmit of \6o\ « 1, equation (15) reduces to the linear growth relation discussed in, for example, Carroll, Press 
& Turner (1992; see also Lahav et al. 1991; Lightman & Schechter 1990), 
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-1 



V 2 

Figure 1 illustrates the different growth rates and expansion factors for the void region and background cosmology. 
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Fig. 1. — The growth of the scale-size and perturbations within a i5o = —0.9 void. This assumes a background cosmology of CIm = 0.3, 
Ha = 0.7, and h = 0.7. Note that at early times the growth of structure within the void and in the background universe axe identical but 
growth in the void is suppressed after t ~ 2Gy (z ~ 5). 
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However, voids may reach well into the 6y{tQ) ~ —0.9 regime, thus the linear approximation is not valid and the void 
universe must be treated numerically. To accurately specify the effective cosmological parameters for large voids, we 
numerically integrate the Friedmann equation for a range of values of rj. This integration yields the parameter a and, 
from equation(15), the present-epoch void density contrast 6y{to). By varying r] until the desired underdensity is found, 
we can compute the effective parameters of the void cosmology. 

In Figure 2, we show the relationships between the density contrast of voids and the corresponding cosmological 
parameters which would be needed to model a universe with the appropriate age and dynamics. Deeper voids expand 
faster, thus they have larger Hubble constant Hy. The matter density parameter Om is smaller, both because of the 
smaller average density and because of the larger Hubble constant (recall the definition fl = {8ttGp)/{3H^)). While the 
cosmological constant A is identical in both the void and background, the corresponding density parameter is smaller 
in the void because the void expands faster. Because voids act like very open universes, the latter reflects the phenomenon 
that dark energy is relatively less important in deep voids, which become curvature-dominated quite early. Even with 
a large cosmological constant, structure formation in voids freezes out at early time in similar fashion to a low-density 
A = universe. 

These results show that a simulation run with hy — 0.84, Qy.M = 0.021, and ily^\ = 0.48 would dynamically evolve just 
as the interior of a spherically-symmetric void with 6 = —0.9 that resides within an = 0.3 flat universe. Obviously, the 
density parameters for the different components of the non-relativistic matter {^oMj^baryon) should be scaled relative 
to those of the background cosmology in the same way. Any patch within such a simulation would mimic the interior of 
a void, up to long-range tidal effects. However, as we discuss in the following section, we must also alter the input power 
spectrum of fluctuations. 

3. THE PRIMORDIAL POWER SPECTRUM IN THE VOIDS 

To correctly specify the input power spectrum for a simulation that seeks to model the interior of large voids, we must 
make allowance for two effects: First, the amplitude of the power spectrum of the density contrast must be adjusted 



Background: Background: 
n„=0.3, 0^=0.7, h=0.7 n„=1.0, n^=0.0, h=0.7 




Fig. 2. — The corresponding cosmological parameters for the interior of voids with varying overdensities, as described in equations (15-17). 
All background cosmological models are flat universes. The parameter r/ is the ratio of the void density contrast to the expansion factor at 
fixed epoch, rj = 5v/a. The efTective Hubble constant is larger in voids because they expand relative to the background, while the density 
parameters Hm and Ha are both smaller (if non-zero). 
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relative to the power spectrum for the background cosmology, because the average matter density is smaller in voids. 
Second, because we wish to simulate voids with fixed density contrast, we should remove power on scales larger than the 
characteristic scale of the voids themselves. Here we discuss these constraints on the initial power spectrum in running 
a simulation of isolated void regions. Again, the goal is to set up the initial conditions so that any region within the 
simulation looks like the inside of a large void, except for possible tidal effects near the edges of the void. 

3.1. Small Scales 

Let us begin by discussing the power spectrum on scales smaller than the void. Consider that at the initial time, U, 
the density fluctuations within some region can be given by 



5{r) = j d^kS^e''^-'' , (19) 



where Sj: denotes the Fourier transform of the density fluctuation field. If the initial conditions are Gaussian, as in the 
current standard models for structure formation, then each of the Fourier components is statistically independent. Thus, 
the contribution to the real-space density field from a particular component is 

dp{f) oc pdj^e'^-'' . (20) 

On small scales, linear independence means that the contribution to the density field from small scale modes is not 
dependent on whether we are looking at a region which will evolve into a void or a cluster. In void simulations of the 
type described above, p in equation (20) is the mean density in the void, not the background universe. Thus, we need to 
apply a correction to the amplitude of each Fourier mode, 

«r'(f.) 

where a superscript, (BG), denotes a random density component in the background cosmology. Recall that rj = 5^{ti)/ai, 
thus this correction is simply the ratio of the background density to the average density within the void. This correction 
yields the same mass perturbations caused by small wavelength fluctuations as in the background. Note that in the limit 
of very high initial redshift this correction almost completely disappears. 

3.2. Large Scales 

We also need to alter the input power spectrum on scales comparable to and larger than the void scale, to tailor our 
simulations to mimic only the interior of voids with specified density contrast. Observed large voids have density contrast 
5 ~ —0.9 and typical radii R ^ lQh~^ Mpc. If we consider that the entire simulation volume is meant to be a "void 
region" then no power should exist on scales much larger than R. 

The variance of fluctuations on a a scale R is 

a''{R,ti) = j 4.nk''dkP^^^\k,ti)\W{kR)\^ (22) 

where W{kR) is the Fourier transform of the spatial window Wb.{x) over which the fluctuations are sampled (e.g., a 
spherical tophat or Gaussian window) and P^^^\k) is the background power spectrum of fluctuations, P{k)^^'^^ = 
(|(5(A;)p). The effect of modifying the amplitude of the Fourier components is to modify the power spectrum as 

P(/e)=/(fc)p(^«)(fc) (23) 

where f{k) is the squared modulus of the factor that multiplies each Fourier component S{k). On scales comparable to 
and larger than the void, we want to smooth away the very fluctuations that led to the creation of a void with density 
contrast Si, such that the amplitude of suppressed power is 

^UpiR^ti) = J AMk[l - /(fc)]p(^«)(fc,t,)|W-(fci?)|2 ~ (24) 

This suppresses suflicient power that the rms of the suppressed structure is similar to the mean underdensity of the void 
region. 

For example, applying a Gaussian cutoff to large-scale power implies 

/(fc) = l-exp[-(fcAo)'] (25) 

as the smoothing function (where ko is the smoothing scale, selected to satisfy equation (24). The result will be little 
structure in the simulations at wavenumber k < ko- Since we do not wish to run a serious of identically periodic voids, this 
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implies that we should limit the size of the simulation box to L ^ 27r/fco. Of course, limiting L effectively removes power 
on large scales. Note that the cutoff scale fcp corresponds to a smaller comoving scale than the current void diameter, 
because all the wave modes within the void grow with time relative to the background. 

Prom above, the smoothing function should approach f{k) = 1/(1 — aiV)^ on small scales. Together, the large and 
small-scale constraints yield the modified initial power spectrum for our void simulations: 



1 - exp[-(fc/fco) 



211 

(26) 



(1 - a,; 77) 2 

4. A RECIPE FOR COOKING VOIDS 

Suppose we wish to simulate a void with density contrast 6y{to) (which specifies the parameter a) and comoving diameter 
at 2; = of io (in units of h~^Mpc, where h is the dimensionless Hubble parameter of the background cosmology). There 
are three separate steps that require careful accounting for the units of distance and scale factors in the void and background 
cosmologies: (1) generating initial density and velocity fields; (2) setting parameters of the simulation code to evolve those 
fields; and (3) translating void simulation results back into comoving coordinates of the background cosmology (i.e., back 
into real comoving coordinates). Here we review these steps. 

To generate the initial density and velocity fields, note that the fluctuations in the void region are those of a patch that 
has comoving size Lo/a. This is smaller than the z = comoving size of the void region because the void has expanded 
relative to the background, thereby stretching all the wavemodes within it. Thus, the amplitude of the fundamental mode 
of the simulation cube should be set by P{k), where k = 27rQ;/Lo- The initial power spectrum must be generated using 
the background cosmological parameters, with amplitude correct for the specified initial redshift Zi, but with shape and 
amplitude modified by the factor f{k) as discussed above. Scaling the velocities, both in the initial conditions and in the 
outputs is very code-specific, as it depends on the units. 

The input parameters for running the simulation include the void cosmological parameters, as well as the comoving 
box size, initial redshift, and final redshift. If specified in units of h~^Mpc, the comoving box size should be set to 
Ly = Lo{h-u/h). The initial redshift of the void simulation, z,,.,;, must be larger than the true Zi, by (1 + Zy^i) /(I + z;) = a. 

To correctly interpret the results of the simulation, we must account for the difference in scale factors and Hubble 
constants in the background cosmology and the void simulation. This is because, although the void expands in comoving 
coordinates when embedded in the background, it does not, of course, grow in comoving coordinates when it sits inside 
the effective void cosmology. At z = 0, scales in the void simulation are exact when expressed in Mpc (not in h~^Mpc). 
Therefore, to get scales in the typical units of h~^Mpc, we must multiply by the ratio of Hubble constants, L = Ly{h/hy). 
At any other epoch, we must also account for the different growth factors, thus the general expression is 

L = Ly^^ (27) 
hy a 

(recall that, by definition, a^(to) = «)• At the initial redshift, z = Zi, L = L{h/hy)/a, as discussed above. 

5. A MASS-FUNCTION SANITY CHECK 

In developing the formalism for a void-only simulation, we have made a number of simplifying assumptions, and 
objections may be raised as to whether the assumptions of spherical symmetry (demanded by Birkhoff's theorem) or fixed 
total mass (demanded by running the void in isolation) are well-founded. 

As a "sanity check," wc have run several relatively modest 128^ particle simulations of a 5y = —0.9 void region. We 
used Enzo (Norman & Bryan 1999), an adaptive mesh refinement (AMR) hydrodynamics code that includes a grid that 
can adapt automatically to provide high-resolution in regions of interest. Refinement was set to allow resolution down to 
1/512 the simulation box-size (25 Mpc), approximately 50 kpc, and only gravitational forces were used. Wc then 
used the HOP (Eisenstein & Hut 1998) group-finder to identify galaxies, and computed the mass-function of the dark 
matter halos in the void region. The results of this simulation can be found in Figure 3. 

We compare the void-only simulation mass function with the mass function found by Gottlober et al. (2003), who 
simulated a large region of the universe, excised the void galaxies and computed their mass function. For consistency, 
the results cited arc based on their Sy < —0.9 sample. Wc also compute a theoretical mass function based on Prcss- 
Schechter (1974) analysis of the spherical collapse of halos. The theoretical curves plotted represent voids of 6y = —0.5 
(top), and Sy = —0.9 (bottom) using the analytic form of Goldberg et al. (2004). The form of Sheth & Tormcn (2002) 
produces a nearly identical curve, but differs in assumption, since the Goldberg et al. result assumes the fully nonlinear 
(semi-analytic) density contrast of the void, derived herein, and uses the size of the void as an additional Bayesian prior. 

It is clear that at most masses, the void-only simulation produces an excellent fit with both larger simulations and 
purely analytic results. The slight discrepancy at the low-mass end can be attributed to limited force resolution in the 
present simulations. As seen in Fig. 3 a simulation with no force refinement produces a significantly worse fit of the 
low-mass mass function. 

6. DISCUSSION AND FUTURE PROSPECTS 

To summarize, the following are the steps to run a void simulation using our new approach: (1) Select the present-epoch 
underdensity Sy{to) of the voids of interest. (2) Select the parameters of the background cosmology in which the voids 
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Fig. 3. — The dark matter halo mass function of a 128"^ particle, Sv = —0.9 void-only simulation, as described in the text (bold solid). The 
dashed line shows the mass function computed in the void regions of the Gottlober ct al. (2003) simulations, while the dotted lines show the 
theoretical mass functions of Sy = —0.5 (top), and 5v = —0.9 void regions using Press-Schechter analysis. The thin solid line represents the 
same simulation with no force refinement, which produces a notably worse fit to both other simulations, and to theory. 
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reside. (3) Compute the effective cosmological parameters (Ou_m, fly, a, Hy) of the interior of the void, foUowing Section 
2 (this requires numerical integration of the Friedmann equation of the void to compute the appropriate vahie of a). 
(4) Compute the modified initial power spectrum of fluctuations, following Section 3. (5) Run a simulation using these 
parameters and power spectrum, in a box with comoving size L, such that the box size is not a large multiple of the void 
diameter. (6) Rescale the simulation results to adjust for the difference in comoving scales in the voids-only simulation 
and the background cosmology. 

Using only a few extra steps in preparing the simulation, we significantly improve the efficiency of our simulations. For 
fixed computational resources, we gain at least an order of magnitude in mass resolution: The voids of interest fill at least 
~ 40% of the volume of the universe but have ~ 10% the mean density of galaxies, thus a voids-only simulation can use 
the same number of particles to simulate the formation of a factor of ^ 25 fewer galaxies at a fixed level of force refinement 
(because we simulate only the void galaxies of interest). High-level adaptive simulations of larger regions will, at best, 
asymptotically approach the efficiency of the proposed method. Simulating only the void region also affords a small 
increase in spatial resolution. This approach is in contrast to running a large simulation, then excising the underdense 
regions for closer inspection. In the latter, huge computational resources are expended on regions which are uninteresting 
for the analysis of voids, i.e., clusters, which soak up lots of computational cycles. 

This approach for running specialized void simulations is valid for studying the interior of large voids. We intentionally 
smooth out structure on the largest scales. Thus, any contribution of tidal structure to the internal dynamics of voids 
cannot be studied by the current approach. Nor, for that matter, can the large scale internal structure of voids be 
examined, as the edges are likely to be strongly affected by nearby structure. 

The next step in this analysis is to apply this approach to run large hydrodynamic simulations of void regions. We 
will use this technique to make projections of the mass spectrum of galaxies in voids (Goldberg et al. 2004), as well as 
their photometric, morphological, and spectroscopic properties. These simulations will be used to interpret results from 
analysis of large samples of void galaxies being identified from the Sloan Digital Sky Survey (Rojas et al. 2003a, b; Hoyle 
et al. 2003), as well as to guide analysis of future deeper surveys. 

DMG acknowledges support from NSF grant AST-0205080. MSV acknowledges support from NSF grant AST-0071201 
and a grant from the John Templeton Foundation. We thank Fiona Hoyle and Henry Winterbottom for useful conversa- 
tions, and Greg Bryan for use of the Enzo code and significant guidance in its use. We thank the anonymous referee for 
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